Flower-like meristem conditions and spatial constraints shape architecture of floral pseudanthia in Apioideae

Background Pseudanthia are multiflowered units that resemble single flowers, frequently by association with pseudocorollas formed by enlarged peripheral florets (ray flowers). Such resemblance is not only superficial, because numerous pseudanthia originate from peculiar reproductive meristems with flower-like characteristics, i.e. floral unit meristems (FUMs). Complex FUM-derived pseudanthia with ray flowers are especially common in Apiaceae, but our knowledge about their patterning is limited. In this paper, we aimed to investigate both the genetic and morphological basis of their development. Results We analysed umbel morphogenesis with SEM in six species representing four clades of Apiaceae subfamily Apioideae with independently acquired floral pseudanthia. Additionally, using in situ hybridization, we investigated expression patterns of LEAFY (LFY), UNUSUAL FLORAL ORGANS (UFO), and CYCLOIDEA (CYC) during umbel development in carrot (Daucus carota subsp. carota). Here, we show that initial differences in size and shape of umbel meristems influence the position of ray flower formation, whereas an interplay between peripheral promotion and spatial constraints in umbellet meristems take part in the establishment of specific patterns of zygomorphy in ray flowers of Apiaceae. This space-dependent patterning results from flower-like morphogenetic traits of the umbel which are also visible at the molecular level. Transcripts of DcLFY are uniformly distributed in the incipient umbel, umbellet and flower meristems, while DcCYC shows divergent expression in central and peripheral florets. Conclusions Our results indicate that umbels develop from determinate reproductive meristems with flower-like characteristics, which supports their recognition as floral units. The great architectural diversity and complexity of pseudanthia in Apiaceae can be explained by the unique conditions of FUMs—an interplay between expression of regulatory genes, specific spatio-temporal ontogenetic constraints and morphogenetic gradients arising during expansion and repetitive fractionation. Alongside Asteraceae, umbellifers constitute an interesting model for investigation of patterning in complex pseudanthia. Supplementary Information The online version contains supplementary material available at 10.1186/s13227-022-00204-6.

visual cue for pollinators can be achieved in various ways, one of them being formation of pseudanthia (aggregates of inconspicuous florets) that resemble and function like a single flower. Pseudanthia evolved independently among many lineages of angiosperms with diverse morphological patterns reflecting their varied phylogenetic descents [3]. Floral pseudanthia with the highest similarity to single flowers are characterized by flowerdimorphism. Peripheral flowers are enlarged and usually zygomorphic ('ray flowers'), whereas the inner flowers are small and radial. Such a morphology is best-known from the heteromorphous heads of Asteraceae (such as daisies or sunflowers), but highly convergent units are also widespread among other campanulids [4,5] which suggests the existence of specific developmental preadaptations for floral pseudanthia in this particular lineage of flowering plants.

The diversity of pseudanthia in Apiaceae
One of the notoriously understudied plant groups that repeatedly acquired floral pseudanthia is Apiaceae subfamily Apioideae (Fig. 1). This large, cosmopolitan clade of campanulid asterids, uniting over 3000 species is distributed across all continents except Antarctica. It comprises numerous economically important crops and aromatic herbs, such as carrot (Daucus carota), parsley (Petroselinum crispum), celery (Apium graveolens), cumin (Cuminum cyminum) and coriander (Coriandrum sativum). The diversity of pseudanthial forms in apioids results from fractal-like architecture of their unique complex umbels, i.e. umbels producing small umbels (umbellets). Pseudocorollas in Apioideae can develop around the entire umbel (umbel-centred promotion) forming a single pseudanthium (Fig. 1B, E), or around each umbellet (umbellet-centred promotion) resulting in several pseudanthia (Fig. 1A, C and F). Ray flowers of umbellifers are also diverse with respect to their pattern of zygomorphy that can encompass different numbers (and parts) of petals. These are often deeply winged due to apical inflexion of their primordia (forming a lobulum inflexum [6,7]). The petal is symmetrical when the wings are equally developed or asymmetrical when the development of one wing is inhibited. Froebe [8] described three different types of zygomorphic pattern formations. The Orlayatype has one enlarged, symmetrically winged petal (Fig. 1C). In the Artedia-type, two asymmetrical petals are enlarged and mirror images to each other (Fig. 1D, E). The Coriandrum type combines both forms by integrating three petals, a symmetrical one in the middle and two asymmetrical petals on its sides (Fig. 1A, B). Another feature of Apiaceae is the diverse initiation sequence of floral organ primordia during flower development. In different species, the sequence can be either centripetal, centrifugal with stamen dominance (they appear as first organs) or sectoral with groups of sepal, stamen and petal primordia appearing consecutively [9]. Despite all those information, until now, the developmental patterning of zygomorphy in ray flowers of Apioideae remains undescribed.

Pseudanthia as floral units
The evolution of biological complexity is frequently based on build-up of simple structures into iterative, compound arrangements [10,11]. This process is apparent in flowerbearing shoots which remarkable diversity is based upon a well-structured, modular architecture. In most angiosperms, an individual flower serves for the basic architectural module that can be born directly on a reproductive shoot (solitary flower) or as a part of a branched system, called inflorescence. An inflorescence meristem (IM) arises from a shoot apical meristem (SAM) and shares many of its qualities, such as apical growth resulting from the maintenance of stem-cell activity (caused by the maintenance of WUS-CLV3 regulatory loop) and ability to acropetally produce axillary meristems via segregation [12]. The self-perpetuation, growth and branching of IM ends with formation of ontogenetically determinate floral meristems (FMs). Depending on the timing and position of FM initiation, inflorescences can develop into simple determinate (botryoids) or indeterminate (racemes) units or their compound equivalents-botryoids of botryoids (compound botryoids/panicles) or botryoids and racemes of racemes (compound racemes).
Traditionally, all pseudanthia were termed inflorescences, as they are composed of numerous florets. However, recently, it has been recognized that alongside inflorescences (originating from IMs), individual flowers might develop as part of various lineage-specific, fundamental modules that repeat themselves on a reproductive shoot [11]. Some of these modules, collectively referred to as floral units, show flower-like developmental qualities. Their ontogenetically determinate meristems (floral unit meristems, FUMs) lack apical activity and instead expand, creating space for centripetal or centrifugal subdivision of submeristems [12,13]. This process, known as fractionation, is controlled by local auxin maxima and proceeds until the entire surface of the FUM is used which implies that the initial size of the meristem and its intrinsic spatio-temporal constraints play substantial role in the patterning of floral units [14,15].
FUMs are frequently associated with floral pseudanthia and constitute basic reproductive modules in several campanulid lineages, including umbels of Apiaceae [12], heads in Asteraceae [16] and Caprifoliaceae subfamily Dipsacoideae [17]. Unfortunately, our knowledge about the genetic patterning of floral pseudanthia is scarce and Fig. 1 Simplified cladogram illustrating phylogenetic relationships of Apiaceae subfamily Apioideae and independent instances of origins of floral pseudanthia (red, modified after 58). Illustrations on the right visualize the diversity of species sampled for morphological investigation (with arrows indicating the clade in which the species is classified), including patterns of promotion and types of ray flower zygomorphy. A Coriandrum sativum with Coriandrum-type symmetry. B Echinophora trichophylla with Coriandrum-type symmetry. C Tordylium apulum with Orlaya-type symmetry. D Tordylium apulum with Artedia-type symmetry. E Artedia squamata showing Artedia-type symmetry of ray flowers. F Scandix pecten-veneris with Scandix-type symmetry outside few model composites, the available data are restricted to genetic mechanisms underlying the elaboration of bilateral symmetry in ray flowers [18][19][20]. As proven by studies on Gerbera and Helianthus, development of the capitulum is governed by genes normally involved in the morphogenesis of single flowers [16]. For instance, its early developmental stages are characterized by the uniform expression of the conserved floral meristem (FM) regulator LEAFY (LFY) that marks it as a determinate structure [21]. Despite this profound change in the quality of the meristem, orthologues of asteracean UNUSUAL FLORAL ORGANS (UFO) genes retained their conventional function related to FM identity [21,22]. The development of showy ray florets is controlled by orthologues of TCP genes CYCLOIDEA (CYC ), normally implicated in the dorsal identity of monosymmetric flowers [23]. Asteraceae possess multiple paralogues of CYC that neofunctionalized creating different expression patterns in ray and disc flowers [24][25][26][27].
The numerous parallels between floral meristems and capitulum meristems indicate that acquisition of floral units may constitute an important prerequisite for evolution of pseudanthia in some plant lineages. The finding of similar molecular patterning in pseudanthia of plants that independently of Asteraceae acquired FUMs would provide further arguments for the formal recognition of floral units and the 'ontogeny-based' concept of inflorescences [13,28]. Here, we present the results of our study focusing on developmental patterns in floral pseudanthia of Apiaceae. Our first aim was to describe how the processes of expansion and fractionation (idiosyncratic for FUMs and FMs) shape differences in the umbel-or umbellet-centred promotion of floral units and in the zygomorphy of their ray flowers. The analysis was based on six species, covering four phylogenetic lineages that independently evolved floral pseudanthia [29]. Our second aim was to analyse expression of three regulatory genes (orthologues of UFO, LFY, CYC ) during the morphogenesis of umbels in wild carrot (Daucus carota subsp. carota) in order to answer the question whether floral units of Apiaceae show 'flower-like' regulation of development.

Morphogenetic patterns in apioid pseudanthia
All analysed pseudanthia show developmental characteristics of floral units at both hierarchical levels (umbel and umbellet). Their patterning proceeds via expansion and fractionation, instead of apical growth and segregation that are to be expected from inflorescence meristems. Interestingly, Apiaceae with umbellet-centred ray-flower promotion and those with umbel-centred ray-flower promotion do not form two separate groups, but represents a continuum of pseudanthial morphologies with intermediate pattern visible in Tordylium brachytaenium. Ray flowers are initiated first, usually as common primordia with their subtending involucellar bracts. They are always perfect and developmentally accelerated in comparison with the inner radial flowers.

Umbel-centred ray-flower promotion
Echinophora trichophylla (Coriandrum-type zygomorphy) Pseudanthia of E. trichophylla develop on the level of the entire umbel (Figs. 1B and 2). The naked FUM ( Fig. 2A) is large (ca. 500 μm) and centripetally fractionates numerous umbellet meristems that are initiated as common primordia with prominent involucral bracts (Fig. 2B, C). The centre of the FUM remains undifferentiated throughout the entire umbel morphogenesis (Fig. 2C, D and E), ultimately giving rise to a plug-like structure (Fig. 2I). Peripheral umbellets develop rapidly (Fig. 2D), overtopping the inner umbellets. The growth of the peripheral involucellar bracts is also accelerated. The pattern of floral meristem initiation in peripheral umbellets is unique. First, two ray flower/involucellar common primordia are fractionated at the abaxial side in some distance to each other. Then, a radial flower primordium appears at the opposite, adaxial side, leading to a distinctly triangular shape of the umbellet meristem (Fig. 2D). Next, the third ray flower meristem arises in the gap between the first two, also as a common primordium with an involucellar bract (Fig. 2D). Central umbellets show a different, spiral sequence of FM initiation (Fig. 2E). The incipient ray flower meristems are spherical and fractionate primordia in the group-like pattern. The first group arises from a common primordium at the abaxial side of the meristem and quickly divides into the primordia of dorsal and lateral petals, the dorsal sepal and its antesepalous stamen (Fig. 2G). Next, a similar group forms at the adaxial side of the meristem, fractionating the second lateral petal, second dorsal sepal and its antesepalous stamen. At the stage of gynoecium initiation (Fig. 2H) the meristem becomes zygomorphic and abaxial sepals elongate significantly into structures resembling 'rabbit ears' marking the onset of Coriandrum-type zygomorphy.
Artedia squamata (Artedia-type zygomorphy) Pseudanthia of A. squamata show an umbel-centred promotion of ray flowers (Figs. 1E and 3) and are virtually indistinguishable from these of E. trichophylla at the stage of umbellet fractionation (Fig. 3A). Its peripheral umbellets are also developmentally accelerated (Fig. 3I) at the umbel meristem which naked centre develops into brush-like structure. The largest difference between the two species is noticeable during the initiation of ray flowers. Initially, two ray flower meristems fractionate at the abaxial side of umbellets, followed by the radial flower meristem at the adaxial side. At the same time, involucral bracts are pressed against the gap between the first two ray flower meristems (Fig. 3B). The two ray flowers in each peripheral umbellet (Fig. 3C) develop rapidly getting oblique position, while their associated involucellar bracts push involucral elements away from the umbel (Fig. 3D). The initiation of floral organs follows the group-like pattern (Fig. 3F). The onset of the future Artedia-type zygomorphy ( Fig. 3H) is established early, when ray flower meristems begin to press against each other, assuming mirrorimaged pentagonal shape (Fig. 3E, F).

Umbellet-centred ray-flower promotion
Coriandrum sativum (Coriandrum-type zygomorphy) Pseudanthia of coriander show distinct, umbellet-centred pattern of promotion (Figs. 1C and 4). Their incipient FUM is small (ca. 100-150 μm) and becomes quickly occupied by few umbellet meristems which fractionate in a weak spiral pattern without associated involucral bracts (Fig. 4A). The development of umbellets proceeds centripetally with rapidly forming raylets (umbellets stalks) separating them vertically (oldest umbellet occupy highest position, Fig. 4B). At this stage, ray flowers are fractionated from umbellet meristems in a centripetal manner as common primordia with involucellar bracts (Fig. 4C), followed by bractless radial flower meristems (Fig. 4D). Involucels grow unevenly; those at the adaxial side of the umbellets are usually smaller (Fig. 4D, E) or even completely obsolete at maturity. The young ray flower meristems are spherical (Fig. 4F) in shape and fractionate floral organ groups (similar to those described in E. trichophylla) in a spiral sequence (Fig. 4G). After initiation of the first three floral whorls, the abaxial sepals enlarge quickly but remain visibly smaller than the petal primordia (Fig. 4H). Lastly, two carpels develop along the abaxial-adaxial axis of the flower indicating inferior gynoecium formation. At this stage, the Coriandrum-type zygomorphy of ray flowers is already established (Fig. 4I).
Tordylium apulum (Orlaya-type zygomorphy) The development of umbellet-centred pseudanthia in this species (Figs. 1C and 5) is almost identical to Coriandrum apart from the presence of small involucral bracts that Scandix pecten-veneris (Scandix-type zygomorphy) The ontogenetic patterns of pseudanthia in S. pecten-veneris (Figs. 1F and 6) deviate from all aforementioned species. Pseudanthia show extremely pronounced umbellet-centred promotion from the earliest stages of morphogenesis (Fig. 6B). While the first of the initiated umbellets fractionates ray flowers, younger ones are still at the naked phase (Fig. 6C). Elongation of raylets separates each of the umbellets into independent developmental units (Fig. 6E, I). Ray flower meristems develop as common primordia with associated involucels that enlarge rapidly into wide, bifid phyllomes (Fig. 6D). The sequence of floral organ initiation begins with stamens (Fig. 6F), followed by petals ( Fig. 6G), and lastly-carpels (Fig. 6H). Sepals are obsolete throughout development. The zygomorphy of ray flowers is established very late and proceeds via elongation of the entire dorsal petal without forming enlarged petal wings (Fig. 6H).

Intermediate form of ray-flower promotion
Tordylium brachytaenium (Artedia-type zygomorphy) An intermediate form of ray-flower promotion can be found in T. brachytaenium (Figs. 1D and 7). Its FUM is similar to that of the umbellet-promoted pseudanthia, but slightly larger (…µm). As the umbellet meristem size is almost equal, this additional space allows for the development of inner umbellets (Fig. 7A). The result is a superficially umbel-promoted pseudanthium with subperipheral umbellets showing smaller, frequently asymmetrical ray flowers (Fig. 7C, D). As in E. trichophylla, peripheral umbellets of T. brachytaenium pass through a triangular stage of development as result of the similar pattern of floral meristem initiation (Fig. 7B). At the early phases of     morphogenesis, ray flower meristems stretch diagonally and become pentagonal, while floral organs fractionate in the already described group-like pattern (Fig. 7E). When carpels begin to differentiate, the outer sepal enlarges with two outer petals (Fig. 7F: sp1, pt1). Each of these petals forms only one wing next to the neighbour flower establishing Artedia-type pattern with two petals being mirror images to each other.

Phylogeny of CYC -like genes
The phylogenetic inference of CYC -like (Fig. 8) gene subfamily supports (with bootstrap = 80) its subdivision in CYC -like genes from early-diverging eudicots (represented by two accessions from the genome of Aquilegia coerulea) and a large, well-defined clade uniting CYC1, CYC2 and CYC3 gene types from core eudicots. In the latter group, only the monophyly of CYC2 genes is Fig. 8 Maximum-likelihood tree of CYC -like genes based on amino acid sequences of conserved TCP and R domains. Major clades are defined with reference to canonical sequences (boldface). The sequence of DcCYC is marked with red. For visualization the tree was rooted with accessions of CYC-like genes from Aquilegia coerulea (early-diverging eudicot). Bootstrap values < 70% were not plotted reinforced with bootstrap value of 73%. The lack of such support for the remaining clades may result from ambiguous position of Vitis vinifera sequence XP002275255.2 that is resolved in a polytomy with respect to CYC1 and CYC3 genes.
The CYC2 clade includes the canonical sequences of CYC and DICHOTOMA (DICH) genes from Antirrhinum majus, as well as TCP1 gene from Arabidopsis thaliana. The large expansion of this gene lineage is clearly visible in sunflower (the genome of which sports eleven paralogues of CYC ). With respect to Apiaceae, both coriander and carrot possess two homologues of CYC , but based on available data it is impossible to unequivocally assess if they arose as a result of two independent or a single ancestral duplication. The carrot's sequence XP017241142.1 constitutes a sister group (with bootstrap support of 100%) to two identical (with respect to TCP and R domain) paralogues from coriander. The second paralogue form carrot's genome-XP017217919.1-is highly divergent from all umbellifer genes belonging to CYC2 lineage.
The phylogenetic reconstruction of CYC -like genes recovered multiple duplication events in asterid CYC1 lineage. The largest expansion occurred in sunflower (six paralogues), followed by carrot (four paralogues), coriander (three paralogues) and snapdragon (three paralogues). A single duplication is also recovered in black cottonwood, while the remaining rosids (Beta vulgaris and Arabidopsis thaliana) possess only one CYC1 gene. The already mentioned XP002275255.2 accession from Vitis vinifera may also belong to this clade. On the other hand, two paralogues from CYC3 clade were identified in most analysed species (cottonwood, wine, snapdragon, carrot and coriander). In this scenario, the single CYC3 genes in Beta vulgaris and Arabidopsis thaliana would arise because of secondary loss. Again, several events of duplication were recovered in sunflower with seven paralogues found in its genome.

Gene expression patterns
The morphogenesis of umbel in carrot proceeds accordingly to other species with umbel-centred pattern of promotion (A. squamata and E. trichophylla). Its incipient FUM is relatively large and flat (Fig. 9A-1) and peripheral umbellet meristems are developmentally accelerated (Fig. 9A-2). Ray flower primordia are fractionated simultaneously with involucellar bracts at the adaxial side of umbellet meristems and enlarge faster than those of radially symmetrical florets (Fig. 9A-3). The pattern of future Coriandrum-type symmetry becomes apparent at the stage of carpel initiation (Fig. 9A-4).
DcLFY is expressed in the mantle zone of the entire FUM during earliest stages of umbel development (Fig. 9B-1). Afterwards, the gene can be detected in fractionating umbellets and flower meristems, as well as in the primordia of involucellar bracts (Fig. 9B-2). The expression of DcLFY is maintained in floral meristems and when floral organs begin to form (Fig. 9B-3) the gene's transcripts can be localized in sepals and petals but not in stamens and carpels (Fig. 9B-4). Contrary to DcLFY, transcripts of DcUFO are absent from the first ontogenetic phases (Fig. 9B-1, 2, C-2) and can be found only in the common primordia of sepals, petals and stamens and later in sepals primordia (Fig. 9C-4).
Based on the phylogenetic inference, we identified two paralogues nested within CYC2 clade (Fig. 8). In their case, the in situ hybridization experiments were focused on the latest developmental stage in which the asymmetry of ray flowers becomes apparent ( Fig. 10A-1). After all organs are fully formed, DcCYC transcripts are found only in ray flowers ( Fig. 10B-1), especially in the enlarged abaxial sepals and medial part of abaxial and lateral petals (Fig. 10B-2). In radial flowers, a weak expression can be traced to one or two peripheral layers of cells surrounding the gynoecium and medial part of the style (Fig. 10B-3). Although the second paralogous gene-Daucus carota XP017217919.1 (Fig. 8)-was expressed in dissected buds (which allowed for synthesis of probes) it could not be detected in developing umbels during in situ hybridization assays, despite several attempts with multiple variants of the basic protocol.

Discussion
For a long time, internode inhibition in simple inflorescences was considered a universal path leading to the evolution of pseudanthia. For instance, capitula of Asteraceae were usually interpreted as derived from simple [30] or compound racemes [31], umbels [32], or spikes [33]. On the other hand, different forms of thyrses (inflorescences with primary racemose and secondary cymose branching) were suggested as the underlying architecture of the umbels in Apiaceae [34,35]. Only recently, the detailed morphological reinvestigation of meristems and increased knowledge about genetic patterning led to the recognition of floral units [12,13,16]. Such multiflowered structures arise from determinate meristems resembling those of single flowers and, thus, are hardly comparable to conventional model plants which inflorescences can retain indeterminate growth for at least some time [36]. The ontogenetic prolongation of FMs, which we hypothesize as the possible pathway towards FUMs [37], may facilitate the co-option of various conserved developmental regulators (such as CYC/ TB1 genes) from the level of individual flower to the level of entire floral unit and constitute an important preadaptation for the evolution of floral pseudanthia. As floral units are morphogenetically determinate, they cannot continuously segregate new submeristems/primordia due to apical growth. Their patterning is thus highly dependent on the space created by the initial and ongoing expansion of the FUM. Additionally, as apioid pseudanthia develop within the envelope of bracts and vegetative leaves that press them against the stem, their patterning is subjected to significant mechanical constraints.

Floral pseudanthia in Apiaceae: the unique interplay among spatial constraints and morphogenetic gradients
Meristem geometry and mechanical constraints during morphogenesis can vastly alter plant morphology [38], including the number of organs [39], their proportion [40] and arrangement [41,42]. Our understanding on how growth rate differences between cells influence gene regulatory networks during plant development is far from comprehensive and although gene expression can be induced solely by mechanical pressure (without secondary, auxin-mediated response), the mechanisms of such mechano-transduction remain elusive) [45,46].
The importance of ontogenetic collisions for architecture of complex floral units has been explored with procedural modelling [47] and is vastly exemplified by experimental studies on simple heads and syncephalia of Asteraceae. Capitula primordia respond to wounding with a change of their phyllotactic pattern and only recently it has been proven that these alterations result from the disruption of natural auxin gradients [15]. In Gerbera, the release of spatial constraints from floral primordia adjacent to the wounding site allows for their repatterning into ray flowers. The bisection of a sunflower capitulum conducted by Marc & Palmer [45] almost four decades earlier yielded a similar outcome, i.e. the formation of two pseudanthial units with enlarged peripheral florets induced at the place of cut. Further proof for space-dependent patterning of FUM-derived pseudanthia comes from natural anomalies of syncephalus Asteraceae [46]. Contrary to the umbels in Asteraceae, syncephalia show distinctly enlarged ray flowers only at the periphery of the entire multi-headed unit. However, single capitula in Oedera capensis and Dyssodia decipiens, whose primordia were physically separated from their neighbours, are able to develop ray flowers around their entire margin [48].
Based on the aforementioned studies and results of our investigation, we hypothesize that patterns of ray flower formation in pseudanthia of Apiaceae result from the interaction between peripheral promotion and spatial constraints which increase towards the centre of the meristem (Fig. 11). Umbel shaping is based upon  umbel-centred (C) promotion patterns in apioid pseudanthia. Green circles represent umbellet meristems. Red circles represent ray flower meristems. The overlapping areas between adjacent umbellet meristems indicate spatial constraints the relative size among the incipient FUM and umbellet meristems. Large FUMs give rise to umbel-centred units (Figs. 2 and 3). Spatial constraints are imposed by the numerous peripheral umbellets, which likely delay the fractionation of inner umbellet meristems and development of adaxial, ray flowers in the peripheral umbellet meristems (Figs. 2D, 3D and 11C). Thus, from the very beginning, the entire pseudanthium develops as a single entity, divided in a promoted peripheral and a retarded central part. On the contrary, when the FUM gives only rise to few umbellets, their physical separation caused by sequential elongation of stalks (raylets) releases mechanical pressure on the adaxial sides (Figs. 4, 5 and 6). In consequence, though all umbellets originate from the same FUM, each of them develops independently and forms ray flowers around its entire margin (Fig. 11A). The establishment of an intermediate promotion pattern in umbels proceeds almost identically to that of umbel-centred pseudanthia. However, due to the smaller size of the FUM and the corresponding smaller number of umbellets (Fig. 7A), its expansion goes along with a formation of additional space between the peripheral umbellets (Figs. 7C and 11B) allowing for the development of subperipheral units with smaller, weakly zygomorphic ray flowers (Fig. 7D). Interestingly, in Apiaceae, the promotion pattern might change within the individual plant (Additional file 1). The higher-order umbels (those that develop later on the reproductive shoot) are usually smaller and produce fewer umbellets. Species with umbel-centred promotion in strong terminal and firstorder umbels can thus produce units with intermediate promotion in higher order as can be seen in carrot (Additional file 1).
The influence of bracts on the polarity and growth of floral meristems is frequently overlooked in developmental studies [36,49,50]. It is widely acknowledged that size and position of bracts might be a source of spatial constraints that influence the shape of floral meristems and floral organ initiation patterns [51][52][53][54][55][56]. In Apioideae, the presence of bracts is highly variable but all species with floral pseudanthia have well-developed involucels [29]. They originate from FUMs as common primordia with ray florets and are thus not a part of the plant's foliage ( Figs. 2A and 3A). The peripheral promotion stimulus acts on those common primordia and causes their simultaneous outgrowth. In the result, similarly to ray flowers, apioid involucels are subjected to spatial constrains that depend on their position within the umbel and overall promotion pattern. In umbellet-centred pseudanthia, the influence of spatial constraints is noticeable in the developmentally retarded adaxial involucels which develop in direct contact with neighbouring umbellet meristems (Figs. 4, 5 and 6). In carrot, the high-order umbels with intermediately promoted ray flowers show a similar intermediate promotion in bracts (Additional file 1). The involucels of peripheral umbellets are large and distinctly pinnatisect, while those of central umbellets are noticeably smaller and needle-like. In subperipheral umbellets, enlarged bracts occur under weak ray florets and are also asymmetrical and less developed in comparison to those found in peripheral umbellets. Besides being subjects of spatial constraints, bracts in Apioideae can also create mechanical pressure themselves. In Artedia squamata, involucre-derived spatial constraints (Fig. 3B) on peripheral umbellet meristems cause retardation of one of the abaxial FMs, as well as changes in the geometry of its neighbours (Fig. 3C). The mechanical forces acting on the sides of ray flower meristems inhibit the enlargement in one of their lateral and dorsal petals, leading to the establishment of Artedia-type zygomorphy. In Tordylium brachytaenium, a similar effect is achieved by the proximity of adjacent FMs ( Fig. 7C and E) which press against each other. The pressure-dependent shift in ray flower meristem symmetry (from radial to zygomorphic) is also apparent in Echinophora trichophylla (Fig. 2H), however its occurrence at late developmental stages-after petals, sepals and stamens are well-developed-does not cause a distortion in the patterning of Coriandrum-type zygomorphy. This observation implies that ontogenetic spatial constraints may have different effects on ray flower morphogenesis, depending on the timing of their occurrence (early or late shift in symmetry sensu Naghiloo, [57].

Floral unit meristems-an important preadaptation for pseudanthia?
The negative feedback loop of CLAVATA3 (CLV3) and WUSCHEL (WUS) constitutes a key genetic component of stem-cell activity, accounting for self-perpetuation of IM and formation of cellular pool necessary for proper development of flowers [58]. While WUS promotes cytokinin activity in the central zone of IM and incipient FMs, CLV3 restricts meristems' size by preventing the build-up of excess cells [59]. The disruption of WUS-CLV3 loop is necessary for cell differentiation and organogenesis. In Arabidopsis thaliana, when a certain size of the floral meristem is reached, WUS acts with LEAFY (LFY) to activate its own repressor-AGAMOUS (AG) [60]. This process establishes the determinacy of the future flower. Other transcription factors that confer floral fate include UFO which works in combination with LFY to specify petals and stamens by activation of B-class MADS box genes [61,62].
FUMs differ from inflorescence meristems and instead resemble flower meristems. Similar to FMs, FUMs are characterized by the early determinacy and lack of apical growth which in Asteraceae coincides with uniform expression of LFY [21] homologue in the naked (i.e. undifferentiated) head meristem. According to our results, umbels in Apiaceae share this pattern (Fig. 9B-1). In carrot, the transcripts of DcLFY can be detected at the umbel meristem and throughout the process of repeated fractionation in the umbellet meristems and flower primordia ( Fig. 9B-2, 3). This indicates that floral units, such as flowerheads and umbels, ontogenetically resemble 'matryoshka dolls' , i.e. or 'umbels within umbels and that the simple process of internode inhibition cannot sufficiently explain their unique developmental patterning. Following Claßen-Bockhoff and Frankenhäuser [37], we hypothesize that FUMs might arise from FMs due to disruption of WUS-CLV3 signalling pathway, including the loss of size-restricting function of CLV3 and/or lack of negative feedback from other direct/indirect repressors of WUS. Such a change would result in the expansion of the meristem and creation of space for additional fractionation, however, this hypothesis requires further studies to be ultimately confirmed or refuted.
The evolution of morphological novelties is frequently based on already existing gene regulatory networks (GRNs). As components of a single GRN are interconnected, differences in spatio-temporal expression of particular regulatory gene can potentially affect its downstream effectors, allowing for their redeployment in the novel context. The examples of such evolutionary co-options are widespread in both animals [63][64][65] and plants [66][67][68]. As FUMs might evolve through heterochronic changes in FM genetic patterning, the co-option of multiple components of flower-specific GRN can tentatively explain the origin of some novelties associated with floral units, including formation of pseudanthia. In all hitherto published studies, including our own (Fig. 10), expression patterns that lead to the establishment of ray flowers were recovered for genes that normally participate either in the specification of identity or symmetry of floral organs [18,[69][70][71][72]. Most of those transcription factors, especially CYC/TB1 genes are known to have independently undergone major expansion in several plant lineages with FUM-derived floral pseudanthia [73][74][75][76][77]. The expression of CYC -like genes is also documented in shoot-derived pseudocorollas of Actinodium cunninghamii [78]. Interestingly, petaloid bracts found in several plant lineages with floral units, such as Nyssaceae [79], Cornaceae [80], are patterned by the expression of B-class MADS box genes [81,82] which normally specify the identity of true petals. Moreover, in dove tree (Davidia involucrata), the origin of petaloid bracts can be traced to duplication and neofunctionalization of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 [83]-a universal flowering pathway regulator [84]. Bractderived pseudcorollas are also known to have evolved independently in numerous clades of Apioideae [29], and in some earlier-diverging umbellifer subfamilies,i.e. species of Alepidea and Astrantia from Saniculoideae, Pozoa coriacea from Azorelloideae and species of Actintous and Xanthosia from Mackinlayoideae [85,86]. Although developmental data are scarce in Apiales, based on the occurrence of the umbel as the basic architectural model in reproductive shoots [87][88][89], we expect floral units to have evolved before the divergence of the clade uniting Apiaceae, Araliaceae and Myodocarpaceae [90,91]. Later, the whole-genome duplication in the common ancestor of Apiaceae lead to the expansion of TCP gene family [92] which might have allowed for repeated co-option of newly acquired paralogues into the patterning of floral pseudanthia.

Conclusions
The results presented in this study provide further evidence that the umbel in Apiaceae, the eponym of its traditional family name Umbelliferae or 'umbel-bearing plants' , should be interpreted as a floral unit, i.e. a multiflowered iterative structure with flower-like characteristics. Its development illustrates all characteristics of floral unit meristems (FUM) including determinacy resulting from flower-like genetic patterning and morphogenesis driven by the processes of expansion and fractionation (self-organizing space-dependent process). As exemplified by surveys of various plant groups, these specific preconditions may constitute an important factor that drives the convergent evolution of pseudanthia. In future, we should aim to identify floral units across angiosperms to gain insight into the diversity of their patterning mechanisms. This in turn, may facilitate understanding of the evolution of plant reproductive structures and shed light on the origin of various morphological novelties associated with them.

Plant material
Wild carrot (Daucus carota subsp. carota) was chosen as the model for investigating gene expression patterns in pseudanthia of Apiaceae. It is a biennial with relatively large and easy-to-manipulate umbels and one of the two umbellifers (the other being coriander) with available genomic data [93,94]. All plants were sampled from natural populations in Warsaw, Poland (52.212671 N, 20.985909 E) and Mainz, Germany (49.992750 N, 8.244798 E). Dissected buds were preserved in cold FAA (50% ethanol, 5% (v/v) acetic acid and 3.7% (v/v) formaldehyde) and briefly vacuum-infiltrated several times. Subsequently, after dehydration and clearing performed according to the protocol of Karlgren et al. [95], the tissues were embedded in ParaPlast Plus ® blocks and stored at 4 ºC until further use. Additionally, a sample of dissected buds was chilled on ice and immediately used for RNA isolation. Morphological analysis was conducted on terminal and first-order umbels of six representatives of Apiaceae (Fig. 1) encompassing both phylogenetic and architectural diversity of pseudanthia in the family [29]. The list of vouchers is provided in Table 1.

Morphological investigation
After bud dissection (between 10 and 30 buds for each species), the specimens were immediately preserved in 70% ethanol. The probes were dehydrated in an ascending ethanol-acetone series, then were critically pointdried (BAL-TEC CPD030), sputter-coated with gold (BAL-TEC SCD005), and observed under the scanning electron microscope (ESEM XL-30 Philips). All steps were conducted according to the manufacturer's protocols. The macrophotographs of carrot umbellets were taken using digital stereomicroscope Leica AM4815.
LFY and UFO orthologues (subsequently abbreviated as DcLFY and DcUFO for D. carota), are single-copy genes in vast majority of angiosperms [99][100][101][102]. Their retrieval was relatively straightforward with a use of protein BLAST [103]. To further confirm our identification, we conducted additional phylogenetic analyses with design similar to that described for TCP genes (Additional file 2).

RNA extraction and in situ hybridization
Total RNA was extracted from dissected buds of D. carota with TRIzol ® reagent. The cDNA for probes synthesis was reverse transcribed using gene-specific primers with SuperScript IV First-Strand Synthesis System (Ther-moFisher Scientific). Promoter for T7 or SP6 polymerase used for in vitro transcription of probes was introduced by PCR amplification with primers carrying appropriate overhangs in 5' end (in forward primer for sense probe, in reverse for antisense, see Additional file 3). The gelpurified products were used in DIG-labelling reaction (DIG RNA labelling kit, Roche). Pretreatment, prehybridization, hybridization and posthybridization steps were performed according to Karlgren et al. [95] using 8-10 μm thick rotary microtome sections (Reichert OmS serial sections). Stringency washes were modified after